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Abstract —This paper presents a real time control strategy for 
dynamically balancing electric demand and supply at local level, 
in a scenario characterized by a HV/MV substation with the 
presence of renewable energy sources in the form of photovoltaic 
generators and an electric energy storage system. The substation 
is connected to the grid and is powered by an equivalent 
traditional power plant playing the role of the bulk power system. 

A Model Predictive Control based approach is proposed, by which 
the active power setpoints for the traditional power plant and 
the storage are continually updated over the time, depending on 
generation costs, storage’s state of charge, foreseen demand and 
production from renewables. The proposed approach is validated 
on a simulation basis, showing its effectiveness in managing 
fluctuations of network demand and photovoltaic generation in 
test and real conditions. 

Keywords—Energy Storage System; Model Predictive Control; 
Demand Response; Renewable Energy Sources; Smart Grid. 

Nomenclature 
Bus active power demand 

Active power generation from traditional power plant 
P^ Storage active power flow 

pres Aggregated active power of photovoltaic generators 
C Cost function associated to traditional power plant 

X Storage’s state of charge 

0 Sampling time 

T Set of time slots in the control horizon 
I. Introduction 

Unlike other types of energy, for many years it has not 
been possible to conveniently store electricity in significant 
amounts. This is the reason why a fundamental requirement for 
traditional power systems has been the one of operating power 
production so as to instantaneously meet the continuously 
changing load demand for active and reactive power. In this 
regard, adequate spinning reserve has to be maintained and 
controlled at all the times to assure a reliable operation of the 
network |[T|. 

As challenging targets in power production from 
Renewable Energy Sources (RES) have been established at 
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international levels 0, the share of fiuctuating Distributed 
Energy Resources (DER) is rapidly increasing in the advanced 
countries, so that a significant need for unconventional bal¬ 
ancing power and storage capabilities aimed at progressively 
relaxing the above-mentioned requirement has emerged in the 
last years. Consequently, technological advancements in the 
field of electric Energy Storage Systems (ESSs) have been 
obtained, putting the basis for new operational strategies of 
electricity systems at high, medium and low voltage levels. 

Some works are becoming to appear in the relevant lit¬ 
erature regarding this new research topic, dedicated both to 
theoretical, applicative and economic aspects. ESS integration 
at residential level is investigated in in relation to 

the problem of supporting consumers in the participation to 
demand side management programs, also taking advantage of 
ICT functionalities offered by Euture Internet technology |8E 
[T^. Other works consider higher size of ESS for applications 
in microgrids and distribution electricity grids |TT[fT^ . Also, a 
significant attention of academics is going towards the implica¬ 
tions of storage integration at transmission level. Dynamic for¬ 
mulations of the optimal power fiow problem at transmission 
level have been introduced and analyzed, both on a simulation 
basis (see e.g. |T5| ) and from a theoretical perspective (see e.g. 
GD, considering a simplified network constituted by a single 
generator and a single load, and GZ) where a proof of concept 
is provided taking the IEEE 14 bus system as case study). 
Eurthermore, in (T^ the integration of ESS and renewable 
energy sources is jointly investigated and simulated. Academic 
and industrial worlds are currently working on this topic, with 
the aim of showing the technical and economic feasibility of 
storage integration through national and international research 
projects (e.g. E-Cube G3 and GRID4EU pQ| ) and pioneering 
demonstrator in operation (e.g (D)- 

In this paper a Model Predictive Control (MPC) strategy 
is presented as a potential tool for facilitating the integration 
of a medium/large-scale electric ESS working at HV/MV 
substation level, with the objective of mitigating the effect 
of fiuctuations in the aggregated distribution grid demand 
and in the generation from DER. The MPC methodology 
has been proved to be effective in a number of industrial 
applications p2}j26| . Its working logic appears promising also 
in the present context: based on a feedback from the field 
and a continuous re-optimization which relies on a prediction 
model providing the state of charge over a control period 
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ahead in time, it realizes ESS integration through a closed loop 
optimal control able to manage fluctuations in the boundary 
conditions and mitigate the effect of model inaccuracies; as a 
matter of fact, the proposed method overcomes conventional 
ESS activation strategies relying on thresholds-crossing and 
schedules designed in advance. Beyond the mitigation of the 
net peak power flowing in the substation, which is modeled 
through an equivalent traditional power plant, the optimality 
criteria used at each iteration of the algorithm also consider 
the need of minimizing the distance of the storage’s state of 
charge from a reference value (in fact, it is important to prevent 
the storage from being fully charged or discharged during 
normal grid operation, so that it can be available to absorb 
potential exceeding energy from RES and support the grid 
in case a contingency occurs). Simulation results are focused 
on showing the relevance of the proposed dynamic approach, 
making use of both artiflcially generated test signals and real 
aggregated demand and photovoltaic generation profiles. 



Fig. 1. Reference scenario. 


III. Modeling 


The remainder of the paper is organized as follows. Section 
[represents the reference scenario in which the ESS is aimed to 
be integrated. Section [11^ is dedicated to the formalization of 
the control problem. In Section IV the simulation results are 
presented and Anally, in Section fVpthe concluding remarks are 
drawn. 


II. Reference scenario 

The reference scenario is depicted in Eig. The core is 
represented by a HV/MV substation where DER and different 
kinds of consumers are connected through the distribution 
lines managed by the Distribution System Operator (DSO). 
The net power flowing through the substation is supposed to 
be provided by an equivalent traditional power plant playing 
the role of the bulk power system. Without loss of technical 
generality, an electric ESS is connected to the MV busbar 
of the transformer in the substation, and controlled by the 
DSO in order to guarantee that the net power flowing from 
the transmission to the distribution network remains close to 
a temporal power profile agreed in advance and is sufflcienty 
smooth. This objective is fundamental to guarantee a natural 
extension of the traditional optimal power flow task to the case 
of DER integration. As a matter of fact, in the traditional power 
flow problem the bus power demand is considered as a known 
input, due to the fact that purely passive transmission network 
nodes can be associated with predictable consumption patterns. 
Such an assumption is challenged (and hence the objective of 
this work becomes extremely relevant) when a significant share 
of DER is integrated at distribution level, due to the fluctuating 
nature of the resource. 

The proposed scenario is sufficiently general to cover also 
other similar situations, for example the one of a stand-alone 
local power system operated in islanding mode, in which a 
real power station is used to feed a large load (for example 
in a industrial area) with the support of RES and an electric 
ESS. Also, from the implementation point of view, the control 
strategy can be hosted by the SCADA control centre operated 
by the DSO through the communication network provided by 
a telco operator or directly hosted into the ESS retrieving 
relevant field data making use of DSO electric infrastructure 
via Power Line Communications (PEC) pT] . 


The proposed control strategy makes use of the MPC 
concept p^ . A receding control horizon is established. At 
each time period the storage’s state of charge is measured 
and taken as initial condition for a proper dynamic prediction 
model. Eorecasts of the bus demand and aggregated RES 
generation stochastic processes are supposed to be known and 
updated over the control horizon. At each MPC iteration, an 
open loop optimal control problem is solved, which is defined 
by interconnecting the balancing power problems related to all 
the time periods in the control horizon using the prediction 
model. The calculated output is represented by the power 
schedule for the traditional power plant and the storage device. 
Reoptimization is used in order to mitigate the effects of inac¬ 
curate system modeling and take into account the updates in 
bus demand over the time. The control problem is formalized 
in the following (symbols are defined in the nomenclature 
section). 

The following objective function to be minimized is con¬ 
sidered, with the aim of establishing a trade-off between the 
need of minimizing generation costs and that of avoiding large 
excursions in the storage state of charge 

^ = E +m [x{t) - x-f} (1) 

teT 

The first term represents the total generation cost from thetra- 
ditional power plant over the control horizon, where each 
contribution can be differently weighted using the parameter 
a{t). In principle, the cost function is a nonlinear, 

time variant function of the generated power. As customary 
in power systems studies | [25| , p9| , the quadratic function 
= 7 (t)P^^(t) is here considered. In the second 
term of the target function, x{t) — is the excursion of the 
storage’s state of charge with respect to its reference; this term 
is summed up over all the the control horizon using different 
weights /3(t), in order to take bounded the evolution over all 
the considered time period. 

A prediction model is used to compute the storage’s state 
of charge evolution over the control horizon. The considered 
dynamics is 


x{t) = x{t — 1) — 0P^(t) \/t e T 


( 2 ) 












IV. Results 


The control P^{t) can be seen as the setpoint established 
by the proposed controller to drive an inner control loop locally 
working at storage level, equivalently to the scheduled power 
generation setpoint P^{t) for the prime mover governor of the 
traditional power plant. However we remark that differently 
from synchronous machines dynamics, storage dynamics is 
included in the proposed control problem formulation, since 
it is not possible to apply indefinitely the same setpoint to the 
storage without incurring in a saturation. 

Regarding the overall behavior of the bus, the balance be¬ 
tween demand and supply over the control horizon is included 
in the problem formulation through the following algebraic 
equation 

P\t)- P'’®* (t) =P%t) + P%t) \/t€T (3) 


A. Overview 

The specific case study under simulation considers an 
HV/MV substation with two transformers, each one rated at 60 
MVA. The equivalent traditional generator is characterized by 
unlimited nominal active power output. The ESS has 12 MWh 
of battery capacity and 6 MW of nominal power input/output. 
Initial conditions and control system references have been 
chosen as = 0 MWh and = 6 MWh respectively. 
The considered sampling period is 0 = 5 m (minutes), while 
the control horizon has been set to 2 hours, unless differently 
indicated. Parameter 7 characterizing the cost function of the 
power plant has been set to 7 = 1 , while the weights appearing 
in the target function Q have been chosen as independent from 
time as (T = 1 and /3 = 5. 


establishing that the difference between the foreseen bus 
power demand and generation from renewables is covered by 
the traditional power plant and storage. Finally, as customary 
in power flow studies, proper box constraints are introduced 
to keep physical variables within their acceptable operating 
range. In this case, the restrictions to be considered are on the 
active power exchange and the storage’s state of charge, which 
are limited as follows 

P3 < p3{t) <P^ MteT 

P^ <P^{t) <P^ WteT (4) 

X < x{t) < X \/t e T 

Finally the real time balancing problem can be stated as 
follows. 

Problem. For a given power bus with known demand P\t) 
and RES generation patterns P^^^{t), generation cost function 
C{P^{t)p) and reference storage’s state of charge x^^^, solve 

min F (5) 

pg ps 

subject to the storage dynamics © with initial conditions x^, 
and the constraints 0(0 over the control horizon T. 

The problem previously defined can be easily put in the 
form of a Quadratic Programming (QP) problem, which can 
be written as ^ 

min + Fu ( 6 ) 

^ineq^ — ®ineg (^) 

AegU = Beg (8) 


< U < U„ 


^min _ ** _ 


(9) 


where vector u of optimization variables is given by the 
collection of variables P^ and P^ at each time of problem 
definition. In particular, 0 can be derived by re-writing ([B 
in matrix form, after having found the explicit solutions of ( 2 ) 
(i.e. after writing x(t) Vt G T in ([^ as function of the current 
state X and of the control variable P^ up to time t). In the 
same way, ^ can be put in the form of equality constraints 
while the sets of constraints 0 can be put in the form of 
inequality constraints 0 and box constraints 0 - 


Simulations have been performed using a Macbook 5.2, 
Intel Core 2 Duo, 2.13 GHz, 5GB RAM 800 Mhz DDR2 
computer, running Apple OSX Yosemite (v. 10.10.1). The 
control framework has been built in Matlab R2012a 64 bit, 
and the QP problem at the basis of the MPC has been solved 
at each iteration by using the commercial solver IBM IFOG 
CPFEX vl2.6. The optimization problem to be solved at each 
iteration consists of 60 variables. The average detected time 
to complete an iteration has been 0.064 s, including the pre¬ 
processing and post-elaboration phases. 

Three test cases have been simulated for providing a proof 
of concept in sample conditions and using real power profiles, 
dedicated to the assessment of control system performance 
at steady state and in response to fiuctuations in the bus 
demand and power production from DER, also considering the 
duration of control horizon as testing parameter. The third test 
is performed using real values of load power and renewable 
energy generation. 


B. Test case 1: fluctuation of demand 

In the first simulation the foreseen aggregated demand 
is characterized by a base load of 50 MW and a peak is 
introduced by adding a Gaussian function centered at 5 pm, 
having an increasing amplitude from 10 % to 50% of the base 
load and a = 4, by which the deviation of the resulting 
network demand from the base load out of the period 4 — 6 pm 
is less than 1% (Fig|^). DER production is not included in 
this first simulation. Starting from the fully discharged storage 
condition (see Figj^, hour 0 ), the system then evolves towards 
a steady-state in which the storage state of charge is close to the 
reference value of 6 MWh (the state of charge is 5.767 MWh 
at hour 3, see Fig ^) and the network demand is completely 
covered by the traditional power plant (Fig hour 3). Such an 
equilibrium depends on the trade-off established in the target 
function 0 between generation costs and penalties for storage 
state of charge deviation with respect to the reference value. 

As soon as the control horizon passes the beginning of 
the perturbation, the storage begins to accumulate energy, 
increasing the battery’s state of charge (see Fig|^. Then, as 
the maximum of the peak demand enters the control period, 
energy is released by the storage, reaching a maximum in 
correspondence of the maximum of bus demand (Figj^). 








Fig. 2. Simulation 1: (a) bus active demand, (b) active power generated by 
the traditional power plant. 




Fig. 3. Simulation 1: storage (a) active power flow and (b) state of charge. 


After that, the storage state of charge recovers the steady- 
state value, ready for new interventions. It is worth noting the 
difference between the demand that should be covered by the 
traditional power plant in absence of the contribution provided 
by the controlled storage (Fig|^), and its actual power pro¬ 
duction (Fig|^), being the last significantly smoother than the 
demand. 

Using the same parameters, the second part of the first 
simulation shows how the length of the control horizon affects 
the results. In Fig|^ and Fig|^ the three curves illustrate the 
responses to the load profile characterized by the peak of 
20% (Fig|^) in presence of a control horizon lasting one 
single time slot (5 m, pointed line), one hour (dotted line) 
and two hours (solid line). The smaller horizon has, as result, 
the largest distance between the storage state of charge and 
its reference value (Fig|^), which brings to a weaker peak 
load mitigation (Fig|^). Increasing the length of the control 
horizon, the storage accumulates energy in advance in order to 
intervene more effectively in presence of the maximum load 
(cfr. Fig|^). A control horizon greater than the peak size does 
not produce relevant improvements. 




Fig. 4. Simulation 1: (a) bus active demand, (b) active power generated by 
the traditional power plant when varying the control horizon. 



t[h] 



Fig. 5. Simulation 1: storage (a) active power flow and (b) state of charge 
varying the control horizon. 

C. Test case 2: fluctuation of RES generation 

In the second simulation, the reaction of the control system 
in response to a fiuctuation in the power production from DER 
has been tested, keeping constant the bus demand to the base 
value. For this purpose a base power generation of = 5 
MW has been considered, then an additional peak power has 
been introduced increasing the amplitude from 25% to 100% 
of the base value and a = 4, by which the deviation of the 
resulting renewable generation from the base value out of the 
period 4 — 6 pm is negligible (Fig|^). 

As in the previous case, once the control horizon passes 
the beginning of the power perturbation, the storage release 
energy (Fig|^) and then has a maximum absorption at the 
peak of DER generation (Fig|7^), then mitigating the variation 
of traditional plant’s power production (Fig|^). 

The red arrows, similar to the previous test, illustrate the 
direction in which the active power produced by the RES 
increases. 
















































































Fig. 6. Simulation 2: (a) active power generated by DER, (b) active power 
generated by the traditional power plant. 




Fig. 7. Simulation 2: storage (a) active power flow and (b) state of charge. 


D. Test case 3: real demand and RES generation profiles 



Fig. 8. Simulation 3: (a) bus demand, (b) aggregated active power generated 
by photovoltaic generators, (c) active power generated by the traditional power 
plant. 



t[h] 



t [ h ] 


In the third simulation, real bus demand (Figj^, pQ| ) 
and photovoltaic generation (Figj^, | [3T| ) patterns have been 
considered, taking 2015-02-21 as tested day. Figj^and Fig|^ 
show the results of this test-case in presence of the proposed 
controller (solid line) and without it (dotted line). Fig|8j; high- 
lights the ability of the control system to mitigate fluctuations 
in the power production from the traditional power plant, at the 
cost of an initial increase of the traditional power generation 
in order to charge the storage (cfr. Fig|^). Moreover, Fig[^ 
shows that the ESS control system based on MFC, besides 
mitigating fluctuations of the renewable generation and, even¬ 
tually, of the load, also makes smooth the most demanding 
power peaks (see for example between 7-9 pm). The resulting 
power flow of the ESS is shown in Fig 9 1 , where the variation 
is within —2 and 2 MW. Finally, Fig 9 3 shows how the the 
storage’s state of charge remains near to the reference and 
tends to restore it whenever possible. 


Fig. 9. Simulation 3: storage (a) active power flow and (b) state of charge. 

V. Conclusions 

A MFC strategy has been presented as potential tool 
for integration of an electric ESS in distribution electricity 
grids, with the objective of mitigating the effect of fluctuating 
aggregated demand and generation from DER on the HV/MV 
substation’s power flow. A proof of concept has been achieved 
in simplifled simulation scenarios characterized by test signals 
and real demand/generation patterns. Future works will be 
dedicated to further testing activities aimed at assessing the 
robustness of the method with respect to mismatches between 
the forecast and the actual realization of involved stochastic 
processes, for example using learning technique like the one 
applied in p2| in order to better take into account the intermit¬ 
tent profile of renewable energy sources. Moreover, a natural 
































































evolution of the addressed problem is the introduction of power 
flow equations in the problem formalization for the extension 
of the traditional off-line optimal power flow problem to its 
real time dynamic version with the aim of integrating 
traditionally non dispatchable generation in the operation of 
transmission networks. 
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